Prescriptions in Loop Quantum Cosmology: A comparative analysis 
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Various prescriptions proposed in the literature to attain the polymeric quantization of a homoge- 
neous and isotropic flat spacetime coupled to a massless scalar field are carefully analyzed in order 
to discuss their differences. A detailed numerical analysis confirms that, for states which are not 
deep in the quantum realm, the expectation values and dispersions of some natural observables of 
interest in cosmology are qualitatively the same for all the considered prescriptions. On the con- 
trary, the amplitude of the wave functions of those states differs considerably at the bounce epoch 
for these prescriptions. This difference cannot be absorbed by a change of representation. Finally, 
the prescriptions with simpler superselection sectors are clearly more efficient from the numerical 
point of view. 
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I. INTRODUCTION 

Loop Quantum Gravity (LQG) [l|, one of the most 
promising approaches to unify general relativity with 
quantum physics, has attracted a lot of attention in re- 
cent years. In particular, considerable progress has been 
achieved in its application to symmetry reduced models 
for cosmology, a field known as Loop Quantum Cosmol- 
ogy (LQC) Q. In this context, the analysis of the sim- 
plest (isotropic) cosmological systems [1, Q has led to a 
qualitatively new picture of the early universe dynam- 
ics [(|, where the current expanding universe is preceded 
by a -semiclassical 0, Hj]- contracting one. This promis- 
ing viewpoint opens new windows in modern cosmology, 
resulting for example in a drastic increase of the proba- 
bility for inflation |9| and ensuring the geodesic complete- 
ness of the (isotropic) cosmological spacetimes [lOj . The 
original analysis has been rigorously extended to various 
topologies [ll| and matter contents [H, EH > as well as to 



homogeneous but anisotropic cosmologies [131 ] . Further- 
more, recent years have witnessed growing progress in the 
extension of the formalism to inhomogeneous scenarios, 
particularly to Gowdy models (both in vacuo [3] and 
with matter [HI) and to perturbative frameworks [l6| . 
In addition, the generalization of the formalism to vari- 
ous Bianchi type models [13] has provided a viable hope 
for a general singularity resolution through the Belinsky- 
Khalatnikov-Lifshitz mechanism |18| . On the other hand, 
the reformulation of LQG as a deparametrized theory [lj[ 
enables the application of the techniques of LQC (either 
directly or after a suitable generalization) in the context 
of the full theory, allowing one to check, in particular, 
the robustness of the LQC results. 

Despite the rapid advances and successful applications 
of LQC to systems of increasing complexity, many of the 
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basic aspects of the theory in its simplest setting remain 
to be fully understood. One of them is the ambiguity in- 
herent to the construction of the quantum Hamiltonian 
constraint. The choice of factor ordering (and densiti- 
zation) gives rise to various quantization prescriptions. 
While all of them provide the same physical picture -to 
a high precision- for the kind of states that are usually 
considered in cosmology (admitting an epoch when the 
universe is semiclassical), any possible physical or math- 
ematical difference in other regimes must be investigated 
and discussed. 

Here we address this question using as probe model the 
simplest cosmological system with nontrivial evolution, 
i.e., a flat Friedmann- Robertson- Walker (FRW) universe 
coupled to a massless scalar field. We focus our attention 
on four prescriptions, three of which have already been 
studied in the literature. These are the original prescrip- 
tion used in Ref. [j| (known as APS, from the initials of 
its authors), its corresponding simplification put forward 
in Ref. [2Q], which allows one to describe the dynam- 
ics analytically (known as sLQC, which stands for soJv- 
able LQC), and the prescription of Ref. [2lj (denoted 
as MMO, again from the initials of the authors). This 
latter prescription is known to significantly simplify the 
physical Hilbert space structure, asserts rigorously the 
generality of the bounce paradigm, and leads to a unique 
Wheeler-DeWitt (WDW) limit in each of the superse- 
lection sectors that are (anti)symmetric under parity re- 
flection. The fourth prescription that we are going to 
analyze is a simpler version of this MMO prescription. 
In this article, we discuss the analytical and numerical 
implications of the application of each of these prescrip- 
tions to construct the quantum model, and we compare 
the details of the physical picture that they provide, in- 
vestigating them in fully quantum (not sharply peaked) 
states. In particular we will show that, while the expecta- 
tion values of certain natural observables show negligible 
discrepancies, physical differences between prescriptions 
actually do exist, making them detectable, at least in 
principle. This fact has important consequences for any 
effective or semiclassical treatment, because it shows that 



2 



the choice of representation and the details of the quan- 
tization procedure actually can have an imprint on the 
dynamics and have to be taken into account. 

This manuscript is organized as follows. In Sec. |TT] 
we briefly describe the classical system and the quantum 
framework. The prescriptions analyzed in this article and 
their main properties are presented in Sec. IIIII Our nu- 
merical methods and results are explained in Sec. lIVI and 
Sec. El Finally, Sec.lVIIcontains a general discussion and 
the main conclusions. In addition, an appendix present- 
ing the WDW quantum counterpart of the considered 
model is included. 



The only nonvanishing constraint that remains after 
the gauge fixing is the Hamiltonian one: 

C(N)=N(C gr + C <jl ), (2.3a) 

2 

C gr = -^c 2 v^H, C^SirGj^-. (2.3b) 

On shell (C gr + C$ = 0), it completely determines the 
dynamics of the system. 



B. Quantum foundations 



II. CLASSICAL AND QUANTUM 
FRAMEWORK 

Let us first remind the construction and basic proper- 
ties of the model in LQC. The foundations and specifica- 
tions of this model have been discussed in Refs. [a, l22j . 
In particular, the APS prescription is described in Ref. 

Details about the othe r q uantization prescriptions 
can be found in Refs. 0, [2l[. We will briefly review 
them all, focusing on those steps where the prescriptions 
differ. 



A. The classical spacetime 

The flat FRW spacetime admits an orthogonal foliation 
by spatial homogeneous 3-surfaces £t (parametrized by 
t). Its metric can be written in the form 

g = -N 2 (t)dt 2 + a 2 (t)°q, (2.1) 

where N is a lapse function, a is the scale factor, and 
°q is a fiducial Cartesian metric, constant in comoving 
coordinates. 

The canonical description derived from the Einstein- 
Hilbcrt action requires integrating the Lagrangian and 
Hamiltonian density over £ t . To avoid divergences of 
the integrals when E t is noncompact, one introduces an 
infrared regulator restricting the integration to a cubical 
cell V (again constant in comoving coordinates). The 
geometrical degrees of freedom are coordinatized in the 
phase space by the Ashtckar-Barbcro connections and 
triads, which, owing to the isotropy of the system, can 
be gauge fixed to the form 

Ai^cV-^Sl E?=pV- 2 / 3 5?, (2.2) 

where V a is the volume of V with respect to °q. Then, all 
the information about the geometry is captured in the 
canonical pair {c,p} = 8tt"/G/3 (where 7 is the Immirzi 
parameter, fixed as explained in Ref. [23[). The matter 
degrees of freedom are described by the field 4> and its 
canonical momentum p^, such that {0,p</>} = 1. 



In order to quantize the system, we apply the Dirac 
program, first representing the classical degrees of free- 
dom as operators, momentarily ignoring the constraint 
(kinematical level). The physical description is then ob- 
tained by imposing the constraint quantum mechanically. 



1. Kinematics 

The kinematical quantization is performed in two 
steps, each of them with a different approach. For the 
matter content, we apply a standard Schrodinger repre- 
sentation. The matter sector of the kinematical Hilbert 
space is identified with H^ in = L 2 (R, def)), spanned by 
the basis of generalized eigenstates (<j)\ of the field oper- 
ator cj). The elementary operators are (f> (which acts by 
multiplication in this representation) and p^, = —ihd^. 

In turn, the geometry is quantized adopting the meth- 
ods of LQG (see Ref. [22| for details). For an isotropic 
model, the standard holonomy-flux algebra can be re- 
stricted to holonomies along straight edges and fluxes 
across unit squares with respect to °q. As a consequence, 
the configuration algebra Cyl s is an algebra of almost pe- 
riodic functions of c, and the kinematical Hilbert space 
for the gravitational sector becomes "Hf- n = L 2 (R, d/is), 
where R is the Bohr compactification of the real line (with 
Bohr measure d/is). A natural basis for "HFL is formed 
by the eigenfunctions \v) of the triad operator p (which 
can be identified with the flux across a unit square) such 
that p\v) — sgn('y)(27r7£| 1 v / A|w|) 2 / 3 |D), where A is re- 
lated with the spectrum of the LQG area operator. This 
basis is orthonormal: (v\v'} = S v . v '. 

All geometric elements of the system (constraints, ob- 
servables) can be expressed in terms of two types of op- 
erators: i) an oriented physical volume corresponding to 
the cell V, V(V) = sgn(p)|p| 3 / 2 , and ii) the holonomy 
components Af^ = e vc / 2 , for an appropriate choice of 
(i. Actually, owing to the specifics of the quantization 
[5(, the choice that must be adopted in the construction 
of the Hamiltonian constraint is such that fi becomes a 
function of the triad, p(p), and this function is fixed by 
the requirement that the square loop with fiducial length 
jl built by the holonomies (in order to define the cur- 
vature) has the minimum physical area that is allowed, 
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A. This choice corresponds to the so called improved dy- 
namics [5| . The action of these operators on the basis of 



V\v) = sgn^lu^^yAlu), Mn\v) = \v + 1). (2.4) 
Finally, the full kinematical Hilbert space is the tensor 

product Hkin = Wjrfn ® ^kin- 



2. The quantum constraint 

We now express the constraint in terms of the oper- 
ators introduced above. This involves, in particular, an 
approximation to the curvature using holonomies along 
a square loop of physical area equal to A, as we have 
already commented (see Ref. || for the detailed proce- 
dure). As a result, the Hamiltonian constraint (|2.3[) . at 
a lapse function of reference, No, takes the general form 



C(N ) = N Q C gr + B 



,-,2 



(2.5) 



where B is some operator which is diagonal in the basis 
and NoCg r is a selfadjoint, difference operator of 
second order. For both of them, the domain of definition 
is chosen to be Cyi s . 

Depending on the prescription, the operator B may 
involve the inverse volume, which is again quantized via 
Thiemann methods and takes the form: 



sgn(V)\V\ 



\ v 47r7/|, 1 v / A / 

x (^iv-i^-^iv-ii/s^y 



(2.6) 



In order to solve the constraint, it is convenient to 
bring it into the explicitly separable form 



C = 8nG(-h 2 & ® 1 + 1 ® f>%) 



(2.7) 



through a process which is often called change of densi- 
tization |24j |. 

The resulting difference operator O is a nonnegative 
selfadjoint [25J operator on Cyl s . For all the considered 
prescriptions, it is defined on Cyl s as 



e = -7V 2p /(«)7V 2 p -M-2 P f(iW-2ji + fo(v) 
where v\v) = v\v) and 
3?rG. 



f(v) 

fo(v) 



4 
3irG 



{v 2 - 2 - a) + 0{v- 2 ) 
-{v 2 -a) + Q{v' 2 ). 



(2.8) 

(2.9a) 
(2.9b) 



Here, a is a constant whose value depends on the factor 
ordering chosen in each prescription. The action of 
only relates states with support on lattices or semilat- 
tices of step 4, denoted by C e , where e S (0, 4]. Then, we 



can identify sectors T-L E C Hkin preserved under the ac- 
tion of 0, as well as by all the observables considered in 
this manuscript (see the definitions in Sec. Ill B 4| . They 
form superselection sectors. At each of these sectors, the 
quantity e used as a label to characterize them, can be 
thought of as proportional to the minimum physical vol- 
ume. 

By applying the analysis of Ref. (25|, one can show 
that the restriction of to each of the superselection 
sectors has essential and absolutely continuous spectra 
that are both equal to R + . Its degeneracy depends both 
on e and on the prescription, but in all the cases it is 
at most twofold. The spectral decomposition of intro- 
duces a natural basis of generalized eigenfunctions (e||, 
solutions to 

[eel}(v) = 12nGk 2 et(v), (2.10) 

where e|(w) = (e||w). The basis is normalized, 

(e%\e%,)=S(k'-k), (2.11) 

where the delta is defined on a domain R which can be 
either R or M. + , depending on the degeneracy of the spec- 
trum (see the discussion in Sec. IIIip . 

Finally, let us clarify the physical meaning of the sub- 
leading coefficient a in Eq. (|2.9p . One can conveniently 
split the evolution operator into 



= O + se, 



(2.12) 



where SQ is a compact operator that contains all the 
terms 0(v~ 2 ) neglected in Eq. (|2.9p . On the other hand, 
©o has quite a simple closed form when expressed in 
terms of v and the corresponding canonical momentum b 
with {v, b} = 4. Namely, introducing the transformation 



WW) = £ e mb '^{v) 



(2.13) 



and the coordinate change x = ln[tan(6/4)]/2, a simple 
computation yields [25] 



Oo 



-YL-kG 



a + 1 



4 cosh (2a;) 



(2.14) 



The prescription-dependent constant a acquires then 
a neat physical interpretation, related to the potential 
strength (see Ref. [25| for more details). 



3. Physical Hilbert space 

Taking into account that C has an invariant dense do- 
main in each T-L e and is essentially selfadjoint there, one 
can construct the physical Hilbert space T-L p h y , for in- 
stance, by applying group averaging techniques [27l [28j 
to Eq. (|27T) . The result is "H phy = L 2 (R,dk), where 
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the domain R is determined by the degeneracy of the 
spectrum, as explained above. This result is completely 
equivalent to the deparametrization of the constrained 
system by selecting <f> as an internal time That pro- 
cedure provides two sectors, which correspond to positive 
and negative frequencies. Without loss of generality, one 
can restrict all considerations e.g. to the positive fre- 
quency sector. In addition, the deparametrization leads 
to a notion of dynamical evolution (consistent with the 
Schrodinger evolution picture) given by the map 

R 3 0^^(-) = *(•,</>); *G^Sn- (2-15) 
The "time translations" are the unitary transformations 

M- = e^^-<^* », (2.16) 

whose generator is the operator \J\Q\ (defined by the 

spectral decomposition of 0). 

Finally, the lack of symmetry breaking interactions in- 
troduces into the system a large symmetry: reflection of 
the triad orientation, v M> —v. It is then possible to split 
the physical space into two sectors: the symmetric and 
the antisymmetric ones. For the rest of our discussion, 
we choose the symmetric sector. This choice does not 
affect the results. 

The physical states take the form: 



k£R 



dk^(k)el(v)e tu{k)<t> , 



(2.17) 



where cj(k) = VT27rG|fe|, e|(w) are the symmetric ba- 
sis functions (|2.10l) corresponding to the superselection 
sector e, and ^ 6 H p h, y - 

4- Natural observables 

In order to extract the physics and test possible dif- 
ferences between prescriptions, we need to introduce (a 
convenient set of) physical observables on "H p hy The 
unitary mapping (|2.15l) allows us to promote any well 
defined kinematical observable O to a physical one, O^, 
acting on the wave function as follows 



= giVI®IW-W 



OV(v,(/>) 



(2.18) 



For our analysis here, we select a set of observables 
that are frequently used in cosmology: the function of 
the volume ln|«|^ [2{J, the energy density -obtained 
from 



A. 

2V 2 



h 2 


T 




T 






e 




2 


V 


V 



(2.19) 



where the symbol ":" stands for symmetric ordering- 
and the Hubble parameter H^, built from 



H = — 



1 



4i7vA 



(2.20) 



All of them leave the spaces H e invariant. 

These operators will be analyzed in Sec. IIV CI where, 
for certain classes of states, we will investigate the differ- 
ence in their expectation values for the distinct prescrip- 
tions under study. These differences between expectation 
values will be considered significant if they are at least of 
the order of the dispersions of the corresponding observ- 
ables. 



III. THE PRESCRIPTIONS 

Even at the classical level, there exists a freedom in 
defining the densitization of the Hamiltonian constraint 
of the system. This amounts to identifying in Eq. (|2.3a[) 
what is the density weight of the function on phase space 
that provides the constraint and which part simply plays 
the role of a Lagrange multiplier. Equivalently, one can 
define the densitization of the constraint by providing a 
(nonvanishing) lapse of reference No, or at least specify- 
ing its density weight, setting then the constraint equal 
to C(N ), as in Eq. (1231) . One can then regard N/N a 
as the associated Lagrange multiplier. Although these 
considerations have no much relevance on purely classi- 
cal grounds, different choices produce different expres- 
sions after rewriting the constraint in terms of triads and 
holonomies; in particular, there appear slight changes in 
the holonomy dependence owing to the regularization of 
the inverse volume terms (see Sec. IIIB"2j) . 

These differences get even more important at the quan- 
tum level where, in addition to the already mentioned 
effects of "holonomization" , we have the freedom to se- 
lect a particular operator representation for the classical 
constraint, including the choice of order in functions of 
noncommuting elementary operators (commonly called 
the factor ordering ambiguity). 

The above ambiguities have already given rise to sev- 
eral constructions for the geometric operator 0. Here, 
we will focus our attention on four of them, three of 
which have already appeared in the literature. They 
are the Ashtekar-Pawlowski-Singh (APS) prescription 
the solvable LQC (sLQC) prescription [20}, and 
the Martin Benito-Mena Marugdn-Olmedo (MMO) pre- 
scription [2l|. In addition to those three, we will also 
consider a prescription (that we call sMMO) that unifies 
some properties of the sLQC and MMO ones. These four 
prescriptions are described below, case by case, pointing 
out the properties of the evolution operator and the form 
of the physical Hilbert space for each of them. 



A. APS 

This is the original prescription used in the pioneer 
analysis of the (improved) LQC dynamics [f| (see also 
Ref. [3(| for a recent discussion of some aspects of this 
quantization). The choices of density weight and opera- 
tor representation are the following: 
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(i) The densitization agrees with the standard one in 
full LQG. Formally, this amounts to setting N = 1 
in Eq. (|2~5]l . 

(ii) For the constraint, one chooses the algebraically 
simple, symmetric factor ordering: 

r^rK-(N 2f ,-N^)V(N 2 f.-N-2f.)- (3.1) 

(iii) In passing to the form ()2.7[) . one requires that the 
emerging Schrodinger system is strictly equivalent 
to the group averaging of the Hamiltonian con- 
straint in its original form (|2.5|) [3lj |. 

As a consequence of these settings, the operator B is 
given by 

= Y M|l« + i| 1/3 -k-i| 1/3 | 3 , (3-2) 

while, the coefficients of the evolution operator O on the 
symmetric sector of the theory are 

f(v) = [f3(v + 2)}^ 2 f(v)[f3(v - 2)] 1 / 2 , (3. 3a ) 

(3.3b) 

+ (l-5 M )/(v-2)], 

where 

/» = ^|»l||»+l|-|«-l||, (3.4a) 

|| B (, r '. ..^o, 

10, v = 0. 

Similar expressions are obtained in the antisymmetric 
sector (except for the 5 Vi ±4 terms, which do not con- 
tribute in that case). 

Furthermore, the implementation of (fTTT|) imposes on 
the physical wave functions (|2.17[) the constraint 

*(0,<£) = 0. (3.5) 

Therefore, in the precise form presented here for this 
prescription, the state \v = 0) does not decouple from the 
rest of the domain from the start (note that in the orig- 
inal representation of Ref. [|[ the wave function still has 
a nontrivial value at v = [30]). However, the appar- 
ent nondecoupling is only formal inasmuch as the state 
\v = 0) does not contribute to the space of physical states 
anyway. 

The constant a describing the potential term in Eq. 
(|2.14j) takes the value: 

a = a A ps = r- (3.6) 
y 

The structure of the physical Hilbert space depends on 
the superselection sector. We distinguish two cases: 



(i) For e = 0, 2 (in the case e = 0, upon our sym- 
metry/antisymmetry assumption) the subdomains 
v > and v < decouple, and the eigenspaces of 
the (extended) evolution operator are nondegener- 
ate. As a consequence, the physical Hilbert space 
is "Hphy' 2 = L 2 (R + ,dk). This case is called excep- 
tional (following Ref. Q). 

(ii) When e =^ 0, 2, the two triad orientations are in- 
terconnected and the eigenspaces of 9 are two- 
dimensional. The resulting physical Hilbert space 
is "Hp hy = i 2 (R, dfc). We will further refer to this 
case as the generic one. 

This dichotomy affects, in particular, the form of the 
symmetric superselection sectors of the theory. For ex- 
ceptional cases, one can restrict all considerations to 
functions supported on semilattices 

Cf = {±(e + 4n),neN} (3.7) 

and next extend them to fully symmetric domains by 
parity reflection. 

For generic cases the superselection sectors have sup- 
port on full lattices, of the form C e = {e + 4n,n G Z}. 
These lattices are not invariant with respect to reflection, 
and therefore one has to work with the union of two lat- 
tices, C £ U £i- £ , first constructing the state on C e and 
then extending it by parity. This will force us to use dif- 
ferent techniques in Sec. HVAl when identifying the basis 
of symmetric functions e|(i>) numerically. 

One of the unfortunate properties of the generic sectors 
in this prescription is the fact that the physical Hilbert 
space is twice larger than for exceptional sectors. For 
the physically interesting applications (analysis of the 
universes which are semiclassical and expanding at late 
times), one restricts the study in practice to half of the 
physical space, thus spanned by only half of the basis 
functions. As we will see in Sec. IIV A[ one chooses them 
to resemble as close as possible (in certain aspects) the 
basis of the exceptional cases. This is achieved by impos- 
ing additional requirements on the behavior of their ge- 
ometrodynamical limit (the selection procedure and nu- 
merical construction will be presented in the mentioned 
section). The remaining basis elements can be then de- 
fined as the orthogonal completion of the constructed 
subset. 



B. sLQC 

This prescription, first suggested in Ref. [2C| , has been 
proposed to bring the evolution operator into a form as 
simple as possible, so that the study of the system can be 
carried out in a fully analytic way. The density weight 
and operator representation of the constraint are as fol- 
lows: 

(i) The constraint is defined with density weight equal 
to one. Formally, this corresponds to a choice of the 
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type iVo = V/(8ttG), for which the time parameter 
is synchronized with the scalar field. 

(ii) The gravitational part of the constraint is defined 
with the ordering 



(3.8) 



With the above criteria, in particular, the operator B 
in Eq. (|2.5[) is just the identity; so no change of densi- 
tization is needed in the quantum theory to reformulate 
the constraint and attain separation of the geometric and 
matter variables. The resulting coefficients of the evolu- 
tion operator are thus much simpler than for the APS 
prescription and read 

f(v) = ^^T2\\ v \t/\^21 (3.9a) 



fo(v) 



4 

3ttG 



(3.9b) 



is 



which implies that the subleading term in Eq. 

a = a sLQC = 0. (310) 

One of the consequences of expression (|3.9p is the anni- 
hilation of the state \v = 0) by 0. The zero volume state 
therefore decouples from the evolution. One can then su- 
perselect this state by its own and remove the quantum 
counterpart of the singularity from the start. The mech- 
anism of singularity resolution is thus slightly different 
from that of the APS prescription. 

The form of O simplifies also the construction of the 
energy density operator, which in this case is 

ft is then possible to show that the entire spectrum of 
p, and not just its essential part [30| (as in the case of 
other prescriptions), is the interval [0, p c ], where p c is the 
critical energy density @. 

Another convenient property of this prescription is the 
fact that, although in the representation used in Eq. 
(|2. 14|) the operator 9 still includes a nontrivial potential 
(apart from the compact remnant), there exists an equiv- 
alent representation in which it adopts a simple Klein- 
Gordon form 6 = -d\ (see Eqs. (3.16) and (3.17) of 
Ref. HJ). 

The structure of the physical Hilbert space is exactly 
the same as for the APS prescription. Again one has 
two cases: (i) the exceptional case for e — 0, 2, where 
the triad orientations either decouple or the degeneracy 
is removed by the parity symmetry restriction, and (ii) 
the generic case, for all other values of e. The treatment 
is exactly the same as in the APS prescription; in partic- 
ular, the physical Hilbert space is %ph y ' — £ 2 (K + ,d£;) 
in the exceptional cases, and the eigenspaces are non- 
degenerate then, whereas for the generic case one has 
U £ phy = L 2 (R, dk) and a twofold de generacy (although in 
practice we will restrict the study to just half of %p hy , 
given the kind of semiclassical states that we want to 
consider). 



C. MMO 

This prescription was originally motivated by the anal- 
ysis of Bianchi I cosmologies in the LQC scenario mill, 
which give rise to natural proposals that, when applied 
to the isotropic situation, affect the factor ordering in 
the representation of the constraint. The main difference 
with respect to the previous prescriptions is the appear- 
ance in C(Nq) of an operator sgn(u). In the isotropic 
context, the prescription was first introduced and stud- 
ied in Ref. [21], following the analysis of its anisotropic 
counterpart. It is characterized by 

(i) The density weight is the same as in LQG. Formally, 
the lapse can be viewed as Nq = 1. 

(ii) The gravitational part of the constraint is defined as 
the reduction of its analog in the Bianchi I model 
[HI], by identifying the degrees of freedom corre- 
sponding to distinct eigendirections of the metric. 

(iii) This and the presence of sgn(u) allows one to choose 
the following operator representation for the con- 
straint: 



N C„ = - 



A 



Mm - M- 



2p 



sgn(w) : A 



(3.12) 



where A is certain operator that is diagonal in \v) 
and satisfies A\0) = 0. We have used the notation 
[: XY :] = (1/2)[XY + XY]. 

In order to achieve separation of variables, one can 
change the densitization and reformulate the constraint 
by a procedure that is directly inherited from the Bianchi 
I model. Namely, one can first deal with the non-isotropic 
Bianchi I constraint and then reduce the result to the 
isotropic case. In this way, one gets a constraint of the 
form (|2.7p . with coefficients for the evolution operator 
(12. 8p that are given by 

/(«) = ^9{v + 2)g(v - 2)g 2 (v)s + (v) S ^(v), (3.13a) 



fo(v) = ^§g 2 (v){[g(v + 2)s + (v)} 



+ [g(v-2)s.(v)} 2 } ) 



where 



s±(v) = sgn(w ± 2) + sgn(v) 

i | | + I|i _ |i _ Us 

9{v) = 



v^0, 
v = 0. 



(3.13b) 

(3.14a) 
(3.14b) 



The resulting evolution operator has several interesting 
properties: 

(i) The operator is an explicit square, O cx Cl 2 , where 
Cl is a known second-order difference operator (see 
Eq. (7) of Hi|). 
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(ii) The coefficients f(v) and f (v) vanish in the whole 
interval v € [—2,2]. 

The latter implies, in particular, that the states whose 
support corresponds to different orientations of the triad 
are not mixed under the action of the constraint. Fur- 
thermore, one can see that the superselection sectors cor- 
responding to © now have support on semilattices (|3.7[) . 
and the absolutely continuous spectrum of in each of 
them is positive and nondegenerate. Hence, each supers- 
election sector (without any exception) has the structure 
found for the exceptional situations in the other prescrip- 
tions presented above. 

This in turn implies that the exact structure of Hph.y 
not only becomes the same for all the superselection sec- 
tors (in contrast with the previous prescriptions), but 
also that it is technically simpler to deal with it (as we 
will discuss in detail in Sec. HV|) . 

Another interesting property follows directly from the 
nondegeneracy of the spectrum. Namely, all the basis 
elements converge in the limit of large v to WDW exact 
standing waves. This property, which holds always in this 
prescription, is achieved in the previous two prescriptions 
only for the exceptional cases e = 0, 2. For the remaining 
sectors, the discussed limit presents a small (decaying 
exponentially with k) but nonvanishing deviation from 
the standing waves, analogous to the case of a tunneling 
through a potential barrier. 

Finally, the asymptotic expansions of f(v) and f (v) 
for v — y oc give 



«MMO 



D. sMMO 



(3.15) 



In this prescription, which has not been stated explic- 
itly in the literature so far, the density weight and the 
operator representation of the constraint are selected to 
bring together the nice features of the sLQC and MMO 
prescriptions: 

(i) The density weight is chosen as in the sLQC pre- 
scription, so that one directly attains separation 
of the geometric and matter variables in the con- 
straint. 

(ii) The constraint operator is defined as a reduction of 
its Bianchi I counterpart, as in the MMO prescrip- 
tion. 

Again, the presence of sgn(w) and the parallelism with 
Bianchi I allow one to choose an operator ordering such 
that the gravitational part of the constraint is of the form 
specified in Eq. (I3.12j) (but with the operator A differing 
from that of the MMO prescription). With these choices, 



O is given by the following coefficients: 
3ttG 



f(v) = — y/\^\\v\y/[v=2\3 + (v)8-(v), (3.16a) 



fo(v) 



3nG 
16~ 



\v\ [\v + 2|4(u) + \v- 2|s 2 »] , (3.16b) 



with s±(v) defined in Eq. (|3.14ap . The asymptotic be- 
havior of these coefficients leads to a potential term in 
Eq. ([2~T4|l with 



a = a s MMO = 0. 



(3.17) 



This prescription shares all the qualitative properties 
of the MMO one: the spectrum, the decoupling of the 
|« = 0) state, the decoupling of the triad orientations 
for all superselection sectors, as well as all the properties 
following from this (see Sec. IIII Cj) . On the other hand, 
the resulting evolution operator differs from the one of 
the sLQC prescription only by a diagonal operator sup- 
ported on v e (—4,4) \ {0}. For the sector e = that 
difference vanishes. This implies that the system is ex- 
actly solvable for e = 0, while for the remaining sectors 
the difference is just a tiny correction (see the compari- 
son in Sec. [Vj) , which for a noncompact system vanishes 
in the limit in which the regulator is removed. 



E. Measuring the differences between prescriptions 

In principle, the considered prescriptions may lead to 
distinct physical pictures. To analyze this possibility, we 
will investigate in Sec. lIVI the behavior of the set of stan- 
dard cosmological observables introduced above. But, 
independently of the existence of significant differences 
for these observables and the feasibility of its detection, 
it is clear that there are discrepancies in the physical sec- 
tor of the theory and that these cannot be absorbed just 
by a change of representation. In fact, we observe dif- 
ferences both in the exact structure of 'Hphy and in the 
(subleading) potential term of O, characterized by the 
constant a [see Eq. (|2.14l) ]. 

This implies that there indeed exist observables which 
can detect the differences between prescriptions, even 
though they may be not of the greatest interest from 
a physical viewpoint. One observable of this kind is the 
densitized constraint C itself. Actually, given a prescrip- 
tion A corresponding to either APS, sLQC, MMO, or 
sMMO, the constraint Ca will obviously annihilate all 
physical states |*a) € "H^ ly , while it will not do so 
generically for physical states of any other prescription 
B 7^ A (provided the action of Ca can be defined on 
them). Mathematically, this means that 



VA^B, 



V X eH kin : {* A \C A \x)=0, 



(3.18) 



3 X eH kin : (*a|C b |x)^0. 
The difference between two constraints is of the form 
C A - C B = AO as ® 1, AQab = ®a- ®b, (3.19) 
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where AQab is a well defined kinematical observable in 
HE . We can then define a family of physical observables 
A6ab \<t> following the procedure explained in Sec. lIIB4l 
Such family allows us not only to detect the differences 
between prescriptions, but also to pinpoint its variation 
in the evolution. 

To understand the nature of these differences, let us 
note that the operator A&ab can be split like in Eq. 
(|2.12p . which gives in the (momentum) b representation: 

AQab = 3ttG aB ~ aA +5e A ~SQ B . (3.20) 
cosh (2x) 

The compact term ASQab = 5@A — 80b can be ne- 
glected in the limit in which the infrared regulator is 
removed. The only residual difference between the pre- 
scriptions is thus generated by the potential term in Eq. 
(|3.20[) . Owing to the shape of the potential, the max- 
imum difference is expected to occur near the bounce 
point x = 0. Its global effect can be understood phys- 
ically as a slight difference in the dispersion of the free 
Klein-Gordon wave packets 32] by this potential. 



IV. NUMERICAL ANALYSIS 

The different properties of the studied prescriptions 
force us to apply different numerical methods in our anal- 
ysis. The elementary "bricks" from which we construct 
the physical states are (some of) the eigenfunctions of the 
operator 0, which form an orthonormal basis of the grav- 
itational part of the kinematical Hilbert space. We will 
thus start (in Sec. IIV Ap with a detailed explanation of 
the procedure to build this basis. We will then describe 
the procedures to construct the physical states and eval- 
uate the expectation values and the dispersions of the 
observables. This will be done in Sec. IIV Bl and IIV CI 
respectively. Finally, in Sec. IIVDI we will discuss and 
compare the efficiency and precision of the methods used 
in different regimes. The details of most of the numerics 
employed in our analysis can be found in Refs. [il. IH l2lj. 
Here, we will concentrate mainly on aspects that are new 
or have not been sufficiently explored before. Our start- 
ing point is provided by appendices B and A. 2 of Refs. [4] 
and [H, respectively. 



A. Basis construction 

As we have already commented, the operator O has a 
continuous spectrum, and is thus diagonalizable in a ba- 
sis of generalized eigenfunctions [solutions to Eq. (12.10)) 
with an infinite kinematical norm] . Depending on the de- 
generacy of the spectrum, these basis elements are sup- 
ported on semilattices (nondegenerate case) or entire lat- 
tices (degenerate case), and can be determined following 
different procedures. 



1. Nondegenerate eigenfunctions 

This is the simplest situation from a technical point of 
view. It occurs in the exceptional superselection sectors 
(e = 0, 2) of the APS and sLQC prescriptions and in all 
sectors of the MMO and sMMO ones. In all these cases, 
the particular form of the evolution operator O is such 
that the eigenfunctions (e|| (where k £ M + ) are uniquely 
determined by their initial value e|(e). To fix the phase 
of the eigenfunctions, we choose this initial piece of data 
to be positive (see Refs. Hl|). Once the values at 
v > are obtained, the function is extended to v < by 
symmetry requirements: e|(— v) = e k (v). 

Let us now explain how to get ef (i>) for v > 0. One 
can see Q that the asymptotic limit v — > oo of e|(i>) has 
the form 

e%(v) -+ r[e**'- e k (v) + e"** e_ k (v)}, (4.1) 

where <jf k is a phase shift, r is a positive real number, and 
e± k (v) are the generalized eigenfunctions of the WDW 
analog of © (see Appendix \K\ . Given that e± k (v) are 
(delta-)normalized as in Eq. (|A3[) . the normalization 
condition (|2.11[) implies that r = 2 (see appendix A. 2 of 
Ref. Q). The relation between the initial value e|(e) and 
the normalization factor r is a priori unknown and can 
be determined only once the eigenfunction (its v — > oo 
limit) is evaluated. To overcome this problem, we divide 
the evaluation into several steps: 

(i) evaluation of a non- normalized eigenfunction ip^; 

(ii) finding its WDW limit to determine its norm ||i/'fc|| 
relative to the condition (I2.11[) : 

(iii) rescaling the eigenfunction to reach a normalized 
one: ip k (v) n- e%(v) = \\ipk\\~ lr 4>k(v)- 

In the first step, we construct ipk by setting ipki^) = E 
The eigenfunction is then evaluated using Eq. (|2.10|) in 
an iterative process, point by point in the domain Cf HI, 
where / = [e,i>Af]. The boundary Vm ^> k is chosen to 
lie (whenever technically possible) deeply in the regime 
where the corrections to the asymptotic behavior (|4.ip 
are small. In the present simulations, following numerical 
estimations, we fix 

vm PS 4 • min[10 7 , max{100 • k, exp(37r/fc)}] . (4.2) 

This choice ensures also that, for small k, the selected 
interval contains at least one oscillation period. 

The second step is completed by a method analogous 
to the transfer matrix technique used in the proof of Eq. 
(|4.1j) in Ref. 0. Namely, the value of the eigenfunction 
at each pair of consecutive points in its domain of defini- 
tion is represented as a linear combination of the WDW 
basis functions, adopting the form 

^k{v) = J— fk(v) cos[fcx + (f>k{v)], (4.3) 
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where x = lnw, and the v-dependent coefficients fk{v) 
and <fik{v) converge according to Eq. (|4.1[) to their re- 
spective limits f\ and <j>%, respectively, with a rate 



fk{v) 
<l>k{v) 







1 + | 










1 + 







(4.4a) 
(4.4b) 



The normalization is now determined by the identity 
f| = 2||^/)fe||. The limit f| can be evaluated quite easily. 
For that we only need a sequence of points {(v^ 1 ,^)} 
with n £ N, extrapolating numerically the desired limit 
at u -1 — > (see appendix B of Ref. Q)- In practice, 
we choose the sequence of points v n to approximately 
follow the behavior v n ~ 2~ n vo, and use a polynomial 
extrapolation (Neville's method). The specific method 
to evaluate fj? depends on the value of fc, namely: 

1. If kxM > 2ir (large k), with xm = lnwMj the 
WDW limit (being a standing wave) has a wave- 
length small enough as to oscillate at least a few 
times in the chosen domain. Collecting the in- 
formation at the extrema of these oscillations, we 
build a set of pairs {(w" 1 ,^)}. The precise algo- 
rithm to evaluate these pairs is the following: 

(a) we find an extremum of tpk(v), namely, a point 
v n £ C+nl where |^fc(«n)| > |^fc(wn+4)| and 
\ipk(v n )\ > \ipk(v n - 4) |; initially, we look for 
the extremum closest to xm\ 

(b) we extend ipk to the interval [v n — 4, v n +4] via 
a polynomial interpolation of second-order (in 
x); the resulting function has the form (14.3[) 
up to fourth-order corrections; 

(c) given this interpolating function, we deter- 
mine the pair {(w" 1 ,^)} corresponding to its 
extremum; 

(d) we repeat the procedure, searching for the 
next extremum close to the point x n — In 2; 
the procedure is repeated iteratively until we 
obtain a sequence of five points, or we enter 
the region where the loop corrections become 
significant. 

2. If Uxm < 2-7T (small k), the wavelength of the os- 
cillations is larger than xm, and we do not get a 
sufficient number of extrema in the selected do- 
main. Therefore, we modify the procedure ex- 
plained above as follows: for each value of v n , 
with x n = lnw„, instead of searching for an ex- 
tremum, we consider the pair of consecutive points 
(v n , v n +4) and solve algebraically Eq. (|4.3I) to find 
{r k (v),<f> k (v)). 

The procedure used for small values of k is simpler, since 
it does not involve identifying extrema nor interpolating, 
but is less accurate than the one employed for large k's. 



Once the sequence {(v^ 1 ,^)} has been found, the 
limit of f j? is determined using a polynomial extrapo- 
lation (Neville's method) at v^ 1 = 0. 



2. Degenerate eigenfunctions 

This is the generic situation (generic superselection 
sectors) found in the APS and sLQC prescriptions. In 
the basis construction, we follow (with minor improve- 
ments) the procedure presented in Refs. [1, @. As al- 
ready mentioned, the eigenspaces are twofold degenerate, 
but in our analysis we concentrate ourselves on a distin- 
guished one-dimensional family of eigenstates (which can 
be next extended to the full basis via orthogonal comple- 
tion) . The general eigenfunctions are solutions of a gen- 
uine second-order difference equation, and hence require 
the specification of two pieces of initial data, e.g. the 
values at two consecutive points of their support. The 
restriction by parity symmetry does not impose any con- 
straints on the data, since for generic sectors the image 
under parity reflection of a lattice C £ is a different lattice 
£4_ e . This implies in particular that any eigenfunction 
supported on C £ can be extended in a straightforward 
way to C e U C±- e by (anti)symmetry. 

Taking into account all this, we construct the (distin- 
guished half) basis as follows: 

(i) first, we build on C e a pair of auxiliary eigenfunc- 
tions V , fc = ( v ) [again solutions to Eq. (|2 . 10[) ] which 
converge to the WDW basis elements e_i fe i in the 
limits v — > ±oo, respectively; 

(ii) then, after a suitable rotation of their phases in a 
process detailed below, we add these functions and 
(delta-) normalize the outcome. 

In the first step of these computations, we choose the 
domain I = C £ [—vm,vm], where vm is selected as 
in Eq. (I4.2[) . The initial data for ■0 ± are given at the 
lowest (for — ) or greatest (for +) pair of points in /, and 
are set equal to the values of e_^(v) at those points. 
While this construction is not an exact implementation 
of ([!]) above, in practice it approximates it with sufficient 
precision, owing to the quick convergence of the LQC 
eigenfunctions to their WDW limits. 

Once the auxiliary eigenfunctions are evaluated, we de- 
termine their WDW limit at the opposite orientation side 
v — > =poo. Since the initial data are complex, this limit 
does not generally correspond to WDW standing waves, 
and takes the more general form: 



(4.5) 



where af,bf £ R + , whereas af,(3^ £ [0,27r). They are 
all functions of k, but we will ignore this in the notation 
so that it does not get too complicated. The numerical 
analysis shows that the amplitude coefficients and 
grow (approximately) in an exponential way with k. 
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And, on the other hand, the selfadjointness of O implies 
that |a±| 2 - |&±| 2 = 1 

To evaluate these coefficients, we split ip^ 1 into real and 
imaginary parts, denoted from now on by the symbols 
and 9, respectively. Since each of them separately con- 
verges to a standing wave, we can then directly apply the 
technique used in the nondegenerate situation, presented 
in Sec. IIV A II The only complication with respect to 
that case is that, in addition to the norm factors [like 
ft in Eq. (I4.4a|) ]. we also need to find the phase shifts 
<j)%. We do so by constructing the sequences {(v" 1 , 
-analogs of {(v' 1 , r^)}- and by finding their limit when 
v^ 1 — >• 0. As before, the values of 0£ are determined for 
large k by the positions of extrema, while for low k they 
are evaluated algebraically. 

Once we know the limiting coefficients of the four com- 
ponents Sftiy^] and Ss^^], the coefficients af, bf, af, 
and fif- can be easily calculated in terms of them. The 
determined data are then used to construct the desired 
linear combination of the two components ■0 ± . This in- 
volves two aspects: (i) normalization, and (ii) rotation. 

Concerning the normalization, we rescale the function 
using the fact that 4||^ || 2 = \af\ 2 + \bf\ 2 + 1. This 
ensures that each of the two considered components con- 
tributes with the same weight to the final basis element. 
Thus, the final result will have a very similar behavior to 
that of the asymptotically standing waves of Sec. IIV A II 

We then rotate ip^ to compensate for the overall phase 



ported on C £ U £4 



1 



(4.6) 



This stepjiiew with respect to the procedure specified in 
Refs. [J, Q, is convenient to improve the semiclassical- 
ity properties of the physical states constructed with our 
techniques from the final basis elements. 

As a result, we obtain the new two components 

^k( v ) = e ' Xe llV'fe'll -1 ^'^ ( v )- One can see that, in the cor- 
responding limits v — > Too, they behave as in Eq. (|4.1[) 
up to corrections of order (a+) _1 , which is a sufficiently 
good approximation for k ^> 1. 

Finally, these components are added and their sum is 
normalized: 



where 



Zk 



a c e 



1 i>j(v) + 4> k (v) 
V2 y/T+Wk\ ' 



+ a+e 



(4.7) 



(4.8) 



This last quantity comes from the scalar product between 
tpt and •fir, and in the regime k 3> 1 is of the order 
of (a+) _1 . As a consequence, it can be neglected for 
physically interesting states (large k). In such case, one 
again recovers the behavior of e| shown in Eq. (|4.1[) . 

The final step in the basis construction is the sym- 
metrization to get the generalized eigenfunction (e| | sup- 



(e%\v) = -j=[(e%\v) + (e%\-v)}. (4.9) 



B. Physical states: construction and analysis 

In the numerical analysis that was carried out in 
Refs. [I,!!, the physical states consisted in Gaussian dis- 
tributions 



*(fc) = * G (k) 



1 



; -(fe-feo) 2 /(4a 2 ) j fc g 



(2^)V4^ 

(4-10) 

The parameters k and a are related in a simple way with 
the expectation value (p<^) and the dispersion Ap^ of the 
momentum of the scalar field 



(p» = Vl2TrGk , 



Ap^ _ jT 

(P<t>) k 



(4.11) 



In principle, the support of the Gaussian spectral profile 
is the entire real line, thus being directly applicable to 
the cases when the basis elements cover the entire set 
fcel (the degenerate situation described in Sec. IIV A 21 
after including the orthogonal complement of the con- 
structed half basis) . In the cases where the spectrum of 
O is nondegenerate, the Gaussian profiles suffer a modifi- 
cation (owing to the cutoff at k = 0). Therefore, in such 
situations the final shape resembles a true Gaussian only 
for profiles that are sharply peaked, so that kg is large 
compared to a. 

Since we are interested in the study of more general 
physical states, for which the different prescriptions may 
in principle lead to different quantum predictions, we 
introduce more convenient profiles, applicable without 
modifications to both the degenerate and the nondegen- 
erate cases. Specifically, we consider logarithmic normal 
distributions of the type: 



1 



(2 7 r) 1 /4 v ^fc 



-[ln(fc/fe )] 2 /(4a 2 ) 



(4.12) 



with k running over the positive semiaxis. The positive 
parameters ko and a are related now to (p<t>) and Ap^ as 
follows: 



(pV) - Vl2^Gfc e" 2/2 , 



Ap> 

(P4>) 



= \Je° 2 - 1. (4.13) 



We will analyze this two-parameter family of states to 
investigate the discrepancies between prescriptions in the 
regimes where (p^) and Ap^ are of the same order. 

The wave function ^{v, 4>) corresponding to a given 
profile 'l'(fc) can be evaluated directly by performing the 
integral (|2.17p . Obvious technical limitations require us 
to, first, discretize the integral and, second, restrict it to a 
compact domain D in k. For our purposes, it is sufficient 
to choose D as 



[k e- sa ,k e s 



with s e N H 



(4.14) 



11 



As far as s > 7, one can check that the relative error 
in the integration owing to the neglected contribution of 
k £ R\D is less than 10~ 12 . In our simulations, we have 
chosen s = 10. 

For the numerical integration of Eq. (|2.17l) in D, we 
have used Romberg's method (see e.g. Ref. [33j). This 
method is particularly convenient if one wants to restrict 
the number of integrand probing points in D, which is 
the case here, as the evaluation of the basis elements 
is the most numerically expensive step of the process. 
To control the integration precision, we have demanded 
that the difference of the results between the consecutive 
orders I and I + 1 of the polynomial extrapolation (an 
internal component of the Romberg method [HI) satisfy 

ll*r i) -*S ) |lphy<'5||*r i) ||p h y, (4.15) 

where the imposed error bound S varies from 10 -6 to 
1CP 10 , depending on the simulation. To avoid an exces- 
sive cost of time in the integration, we have restricted 
the number of the integrand probing points (forming the 
uniform lattice in D) to 2 12 + 1. 

C. Observables 

We can now proceed to calculate the action of observ- 
ables on the states represented by the wave functions 
constructed in the previous subsection. In particular, 
we can evaluate and compare the expectation values and 
dispersions of those observables. We consider two types 
of observables: those introduced in Sec. Ill B 41 which 
encode standard properties of interest in cosmology - 
namely, ln|{i|^, H^, and p^-, and the observables de- 
fined in Sec. IIIIE1 which are specially suitable to detect 
the differences between the studied prescriptions. 

The dynamical information is extracted by means of 
the Schrodinger picture, where the evolution of a state 
is seen as a mapping between initial data (on a constant 
<f> slice) via the unitary transformation (|2.16[) . In this 
picture, the action of a physical observable is obtained 
from that of its kinematical precursor (see Sec. Ill B 4[) on 
the appropriate initial data slice. 

This fact has been applied in previous numerical stud- 
ies of LQC, starting with Ref. in order to extract 
dynamical data. In our case, in the v representation, 
all the interesting kinematical operators (precursors) are 
either multiplication operators or combinations of mul- 
tiplications and shifts. This simple form allows us to 
evaluate the results of their action straightforwardly by 
making use of the map Vf^(u) i-» hx\v\^^,(v), Eq. (|2.19l) . 
Eq. (|23D|) . and the right formula in Eq. (IXTS)) . 

The expectation values arc then evaluated via the kine- 
matical inner product on H gI , 

(*,|$*)= ]T (4.16) 
vec c nJ 

where, for technical reasons, the summation is restricted 
to the compact region J = [— v m i, v m >] (in the degenerate 



case) or J = [0,u m /] (in the nondegenerate case). In our 
simulations the bound v m > has been selected to lay always 
in the interval [10 4 ,4 • 10 6 ], its specific value varying for 
different simulations. This choice ensures that the error 
caused by the restriction to a compact domain J has a 
subleading contribution compared with other numerical 
errors that arise in the evaluation process. 

To isolate the numerical noise (see Fig.[Tj generated by 
the errors introduced in the evaluation of the basis and in 
the integration of the wave function ^ on each slice, the 
values of ^^(v) entering Eq. (14. 16|) have been modified 
by a filter, namely, whenever |^0(u)| < asupl^/^l the 
value ^4>(v) has been set equal to zero. This prevents 
this type of noise from affecting the computation of the 
expectation values. In our simulations, the value of the 
relative bound a has been selected to vary between 1CP 8 
and 1(T 6 . 

The dispersions have been calculated using the stan- 
dard formula 

(AO} 2 H0 2 )-(0) 2 (4-17) 

for each observable O. 

Finally, since we work in the symmetric sector, we 
note that we can restrict all our considerations to half 
of the support of the wave function; in particular, when 
limited to a compact region, we can restrict ourselves to 
a domain like J, as specified below Eq. (|4.16p . 



D. Efficiency and precision 

Let us now investigate the properties of the numeri- 
cal techniques discussed in the previous subsections from 
the perspective of the computational precision and ef- 
ficiency. The numerical computations necessary to get 
the final results consist of several steps: (i) evaluation of 
the basis elements (e|| (discussed in detail in Sec. II V Al) . 
(ii) integration of the spectral profile to determine the 
wave function for each constant value of <j> (Sec. IIV Bl) . 
and (iii) evaluation of the action of the observables and 
computation of their expectation values and dispersions 
( Sec. IIV C\i . Each of these steps introduces its own source 
of numerical error and presents a different efficiency. 

We start with the first of these steps: the evaluation of 
the basis elements. The comparisons during the simula- 
tions have shown that this step is responsible for most of 
the computational cost, and therefore it is the most criti- 
cal part from the viewpoint of the efficiency. As discussed 
in Sec. IIV Al the actual algorithms and cost depend on 
the degeneracy of the basis, and hence vary significantly 
with the considered superselection sector and quantiza- 
tion prescription. 

In the nondegenerate case (Sec. IIV A 1[) . the calcula- 
tion involves two steps: the determination of the non- 
normalized eigenfunctions ?pk and the computation of 
their norm by finding their WDW limit. The calculation 
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(a) (b) 

FIG. 1. Amplitude ^^(v)! of the physical wave function corresponding to a state with a logarithmic normal distribution. The 
amplitude for different prescriptions is compared both away from the bounce (a) and at the bounce (b). The parameters of the 
profile ^ of this state are fixed by the conditions {p,p} — 100ft and Ap^/{p0) = 0.1. Away from the bounce, the amplitudes 
are indistinguishable up to numerical noise, whereas at the bounce one can observe differences (phase shift) in the interference 
pattern. The noise level clearly depends on the used techniques, something which depends in turn on the degeneracy of the 
spectrum of ©. 



precision depends on the size of the evaluation domain 
chosen for the eigenfunction and on the wave number k 
[see Eq. (|2.17l) ]. In particular, we observe that two effects 
compete: since the eigenfunctions are evaluated via itera- 
tive methods, the evaluation precision decreases with the 
size of the domain, whereas the precision in determining 
the WDW limit increases with it. In that respect, the 
choice of the domain size given by Eq. (|4.2I) provides a 
fairly acceptable balance between these two sources of er- 
ror. It is also worth recalling that, with our conventions 
[e|(e) > 0], there are no ambiguities in the freedom of 
choice for the global phase of the eigenfunctions. 

The degenerate case, as we have seen in Sec. IIV A 21 
is considerably more complicated. First, the procedure 
applied in the nondegenerate case becomes just the first 
step of the evaluation. Even this stage introduces now a 
higher numerical error, because the domain of calculation 
of ipk is now twice larger, and hence the evaluation of the 
eigenfunctions requires twice more iterative steps. Apart 
from that, we observe a significant cost increase since 
we have to evaluate the pair of eigenfunctions ip^ and, 
besides, both %jxr are now complex instead of real. In 
total, the three commented facts amount to an increase 
of 8 times in the computational cost. 

Furthermore, the next step -taking the appropriate lin- 
ear combination of ip^ to form the final basis functions- 
has its own cost (which is linear in the domain size). 
Apart from that, the rotational symmetry of the compo- 
nents tp^ is broken, in the sense that, in order to con- 
struct the appropriate basis vectors, we need to compen- 
sate for the overall phase of the WDW limits of those 
components (14.61) . This step introduces extra complica- 
tions, since the phase itself can be determined only mod- 
ulo 7r. The correct identification of this phase, crucial 



for the subsequent construction of the relevant physical 
states, is nontrivial, and in fact one can check that this 
phase is proportional to A;ln|fc| at its leading order. As 
a consequence, this step in the determination of e|(u) 
introduces an additional numerical error. 

In the integration of the wave function profiles [step 
(ii) above], the use of the high-order Romberg's method 
allows us to restrict the number of evaluated basis ele- 
ments to a manageable amount, while keeping sufficiently 
high numerical precision. The selection of this method 
and of a proper compact integration domain makes also 
possible that both the integration error and the error 
caused by the restriction of the domain can be limited 
to a level where they do not exceed the error generated 
in our previous step (i) of the numerical computation. 
The differences between the degenerate and nondegener- 
ate cases do not require a different treatment. However, 
in practice, the degenerate situation turns out to be ap- 
proximately 3 times more expensive numerically owing 
to two reasons: (a) because the eigenfunctions e\(v) are 
complex in that case, and (b) because the wave function 
has to be calculated for both v > and v < 0. 

The effect and dependence of the overall numerical er- 
ror introduced in the previous steps (i) and (ii) is shown 
in Fig. [TJ For the states analyzed in this article, the er- 
ror stays at the level of 10 -12 in the nondegenerate case. 
The additional complications characteristic of the degen- 
erate case cause the error to grow in those cases by 2 or 3 
orders of magnitude. Nonetheless, all the wave function 
profiles can be integrated with a final relative error which 
does not exceed 10~ 8 . 

The final step (Hi) in the numerical computation in- 
volves algorithms which are common for both the de- 
generate and the nondegenerate cases. However, in the 
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degenerate case, a higher level of numerical noise is visi- 
ble in Fig. [1] This has forced us to conveniently increase 
in this case the value of the relative bound a in the dis- 
crimination filter (see Sec. lIV Cj) . In turn, this happens to 
increase the error in the evaluation of the expectation val- 
ues and dispersions by the same order of magnitude (i.e., 
it increases from approximately 10~ 12 to 10 -9 -10 -8 ). 



V. RESULTS AND DISCUSSION 

We have applied the methods explained in the previous 
section to the numerical analysis of a population of states 
that are given by a normal logarithmic distribution of the 
form (|4.10p . with the value of (p^) ranging from 30ft to 
500ft and the relative dispersion (/S.p<p) / (p^) from 0.05 to 
0.25. The analysis of these states has been carried out in 
the four prescriptions discussed in this article. We have 
analyzed 2 5 different values of e, labeling distinct super- 
selection sectors. The results are displayed in Figs. [TUB] 
At various levels of comparison, we can distinguish the 
following aspects. 

First, a preliminary comparison can be performed at 
the level of the wave function itself. Namely, one can 
compare the probability amplitude -the value \^(v, <f>)\ 
scaled by the square root of the inner product measure on 
"H gr - of the wave function which represents the same state 
[i.e., with the same spectral profile \P(fc)] in the different 
prescriptions. Away from the bounce [see Fig. l(a)| , one 
does not observe any significant difference. However, at 
the bounce [Fig. 1(b) the situation becomes slightly more 
complicated. 



The general shape of the wave function 
(position of the peak, general behavior of the function 
slopes) still does not show any clear distinction; never- 
theless, one observes a phenomenon that actually reveals 
the existence of fine differences. In fact, the interaction 
of the expanding and contracting branches creates an in- 
terference pattern, which can be seen on the downward 
slope of the function in Fig. 1(b) For a chosen superse- 
lection sector and wave functions representing the same 
state, the pattern shows a dependence on the prescrip- 
tion: for different prescriptions the minima and maxima 
of the interference are displaced. Nonetheless, the spe- 
cific shift of the various extrema depends not only on the 
prescription, but also on the shape of the state (spec- 
tral profile), as well as on the superselection sector to 
which it belongs. As a consequence, it cannot be used 
in a systematic obvious way to identify the prescription 
employed, regardless of the state under consideration. 

Another, more physically relevant aspect for compari- 
son has to do with the use of the cosmological observables 
ln|u|0, Hrf,, and p,£. The results are presented in Figs. [2]- 
[5j Analyzing the same state in different prescriptions we 
have found detectable differences between the ex pecta- 
tio n val ues of all the three observables [see Figs. |3(b)| 
SI 5(a)| . These differences are most prominent at the 
bounce and decay quickly as the wave packet enters the 
low energy density regime. For the states investigated 



here, for which (p^) > 30ft and (Ap^,) < 0.25(p^,), the 
differences are nevertheless several orders of magnitude 
smaller than the dispersions of the corresponding observ- 
ables through all the evolution. Those differences depend 
on the degeneration of the spectrum of 0, apart from 
their natural dependence on the observables and the par- 
ticular state under consideration. The situation where 
the highest differences have been observed occurs when 
one compares the expectation values of the energy den- 
sity operator on highly dispersed states with low momen- 
tum (pcj,). Then, the computed differences lay only one 
or two orders of magnitude below the dispersions during 
the whole evolution. In the rest of situations considered 
here, the differences are even smaller when compared to 
the dispersions. 

Among the results presented above, the dispersion of 
the energy density deserves special attention. For all 
the prescriptions, the essential spectrum of this operator 
is the interval [0,p c ] where p c s» 0.81ppi is the so called 
critical energy density. Depending on the prescription, 
the spectrum may also have a discrete part with eigen- 
values exceeding p c , but these play no role in the states 
that represent the cosmological solutions [30]. This fact 
is reflected in the behavior of (Ap^). Namely, for the 
states analyzed in this paper and for the nondegener- 
ate cases (as defined in Sec. lIV A|) . the expectation value 
(ptp) reaches the critical energy density p c at the bounce, 
and its dispersion drops significantly there (in particular, 
it vanishes up to the numerical error for states peaked 
at large p$). In this sense, the states with the spec- 
tral profile (14. 10)) are coherent ones. In the degenerate 
cases the situation is different: we observe that (Ap^) 
decreases near the bounce, but it reaches a positive mini- 
mum significantly larger (at least a few times) than in the 
nondegenerate case. This property clearly distinguishes 
between the APS and sLQC prescriptions, on the one 
hand, and the MMO and sMMO prescriptions, on the 
other hand, at least for the superselection sectors with 
e ^ 0,2. The observed difference might be nonetheless 
related to the particular way of constructing the basis in 
the degenerate cases |34 |. 

The other physical aspect considered in our numer- 
ical analysis concerns the expectation values of the ob- 
servables (AQab) 2 , constructed specially to measure the 
discrepancies between the different prescriptions. The re- 
sults are presented in Fig. [5] As we can see, these expec- 
tation values are nonvanishing. Thus the prescriptions 
are clearly distinct, and the differences between them are 
of course detectable. As we could have guessed, the ex- 
pectation values (and therefore the physical differences 
between the prescriptions) are largest near the bounce. 
Away from it, they decay exponentially. This behavior 
can in fact be proven analytically for all physical states, 
for which (Ap^) is finite, by employing methods similar 
to those applied in Sec. VB of Ref. [§| . Not surprisingly, 
the largest differences are observed between prescriptions 
which lead to a different potential term in the expansion 
(|2.14[) (i.e., to different values of the constant a). An ex- 
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FIG. 2. Dynamical trajectories of H$ (a) and p^, (6), given by the expectation values of these observables on the state of Fig. [T] 
with e — 1 in the APS prescription. 
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FIG. 3. (a) Quantum trajectory of In |?5|^ for the same state and the same prescription as in Fig. [2] (6) Uncertainty in In \v\^, 
for the same state and prescription, compared with the difference between the corresponding expectation values calculated in 
the APS and MMO prescriptions. 



ample of such situation is presented in Fig. 6(a) In the 
case of the sLQC and the sMMO prescriptions [Fig. 6(b)| , 
the difference is many orders of magnitude smaller (in the 
presented case, more than 16 orders), because the poten- 
tial terms in these two prescriptions coincide and the only 
difference is a compact operator, supported only on three 
lattice points near the classical singularity. 



VI. CONCLUSIONS 

In LQC, even in simplest models, the standard ambi- 
guities of the canonical quantization affecting the spec- 
ification of the Hamiltonian constraint and its operator 
representation have led to several quantization prescrip- 
tions. In this paper we have analyzed in detail three of 
those most commonly used in the literature, known as the 



APS @, the sLQC 0, and the MMO [H[ prescriptions. 
In addition, we have introduced a new one, the sMMO 
prescription, which combines useful features of both the 
MMO and the sLQC ones (see Sec. lTlTDi 

Basically, different prescriptions lead to slight differ- 
ences in the evolution operator O that generates the uni- 
tary dynamical evolution in the internal time, whose role 
is played by the massless scalar field. These differences 
make that the use of one or another of the prescriptions 
results to be more convenient in distinct circumstances, 
depending on the particular application under consider- 
ation. 

In particular, the mathematical structure of the physi- 
cal Hilbert space is different for the various prescriptions 
discussed here. In fact, for generic superselection sec- 
tors, the system has a rather more complicated structure 
in the APS and sLQC prescriptions, owing to a twofold 
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FIG. 4. Absolute dispersions (a) and relative dispersions (6) of for the considered state of Fig. [TJ compared with the 
corresponding difference between the expectation values of in the APS and MMO prescriptions. For both relative values, 
one can observe a peak at the bounce owing to the vanishing of (H$); however, the peak in the differences (red curve) is so 
sharp that it is placed between probing points. 




FIG. 5. (a) Relative dispersion of for the same state as in the previous figures, compared with the relative difference between 
the corresponding expectation values in the APS and MMO prescriptions, (b) Comparison between the relative dispersion of 
in the APS and MMO prescriptions, for a "generic" superselection sector. Better coherence properties are observed for the 
state constructed in the MMO prescription. 



degeneracy of the spectrum of 0, whereas in the same 
situations the MMO and sMMO prescriptions (for which 
the spectrum is nondegenerate) provide a much simpler 
Hilbert space of physical states. This fact has a signif- 
icant influence in the efficiency of the numerical tech- 
niques used in the dynamical study of the system, which 
therefore varies considerably from the first to the second 
of these sets of prescriptions. As discussed in Sec. lIVl the 
construction and analysis of the physical states in the de- 
generate cases requires more complicated methods, which 
in turn increase the computational cost and the numer- 
ical error. Although this error is far from critical in the 
computations that we have performed, since the relative 
error grows in the degenerate cases, compared to the non- 



degenerate ones, from approximately 1CP 12 to only 10~ 9 , 
the problem of the time cost is relevant from the numer- 
ical point of view. As discussed in Sec. IIV Al the cost of 
the (most demanding) step, in which the basis of "H p hy is 
constructed, is at least 8 times higher in the degenerate 
cases than in the nondegenerate ones. This shows that, 
whenever the system has to be analyzed numerically, the 
MMO and sMMO prescriptions are much more appro- 
priate. The cost difference becomes particularly critical 
once one tries to analyze more complicated cosmological 
models, like for example Bianchi I |35j |. 

Despite the significant focus on technical aspects, the 
main aim of our investigation has been to identify and 
study possible differences between the considered pre- 
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FIG. 6. Expectation values of the observable (AOab) 2 |</> for A = APS and B = MMO (a); as well as for A — sLQC and 
B = sMMO (6). In both cases, e = 1 and the observable is evaluated on the state with a logarithmic normal distribution 
whose parameters are (prf,) = 100ft and Ap r /,/{p r i,) — 0.1. The state is built in the MMO prescription in case (a), and in the 
sMMO prescription in case (b). The difference reaches a maximum at the bounce and decays exponentially away from it. The 
difference between the sLQC and the sMMO prescriptions is many orders of magnitude smaller than the difference between 
any other pair of prescriptions. 



scriptions on a physical level. To achieve this, we have 
analyzed a two-parameter family of physical states with 
spectral profiles corresponding to a logarithmic normal 
distribution [see Eq. (14. 10)) ] . and without imposing the 
restriction of semiclassicality. For our analysis, we have 
chosen states peaked about low values of the scalar field 
momentum, (p^) < 500S, since the differences are easier 
to unveil in this regime. For these states, we have been 
able to detect differences between the various prescrip- 
tions by observing the interference pattern in the wave 
packet tail at the bounce. The comparison of the states 
built in the different prescriptions, for the same spectral 
profiles, has shown that the patterns are actually shifted 
with respect to each other. This result confirms the exis- 
tence of differences. Nonetheless, it does not allow one to 
straightforwardly deduce which specific prescription has 
been employed, because the commented shift depends 
also on other factors, such as the superselection sector 
and the particular spectral profile of the state. 

In order to address in depth the feasibility of the de- 
tection of differences between prescriptions, further anal- 
ysis has been performed. We have focused it on two 
fronts, discussing the discrepancies in the measurements 
of standard cosmological observables, on the one hand, 
and studying certain quantum operators which are spe- 
cially sensible to a change of prescription, on the other 
hand. 

Concerning the first of these fronts, we have picked 
up three observables of interest in cosmology, namely, 
the logarithmic volume, the Hubble parameter, and the 
scalar field energy density. We have used them to com- 
pare the dynamical (quantum) trajectories of the physi- 
cal states specified above. We have evaluated the differ- 
ences in the expectation values of these observables be- 



tween the different prescriptions, and shown that they are 
several orders of magnitude smaller than the respective 
dispersions. As a consequence, and as far as we restrict 
ourselves to these standard cosmological observables and 
to the considered physical states, the differences between 
prescriptions are not detectable in practice. 

As for the other kind of observables that has been con- 
sidered, we have computed the expectation values of the 
operators (AQab) 2 defined in Eq. (I3.19|) . which essen- 
tially encode the differences between the Hamiltonian 
constraints that correspond to different prescriptions. 
These expectation values are nonvanishing, reach the 
maximum at the bounce, and decay exponentially away 
from it. The nature of the physical differences between 
prescriptions has been well understood (see Sec. IIIIEI) . 
The principal component from which these differences 
arise comes from the potential term of in the volume 
momentum representation. This potential term does not 
coincide for all the studied prescriptions. While the sub- 
leading remnant in O also varies when one changes the 
prescription, this remnant is a compact operator and its 
effect is negligible. This also explains why the smallest 
differences are observed between the sLQC and sMMO 
prescriptions, since the potential term is the same in 
these two cases. 

At this point, it is worth recalling that the studied 
physical trajectories and the measured differences are 
genuinely well defined only if the spatial homogeneous 
slices are compact (in the considered model, of T 3 topol- 
ogy). In noncompact cases, it is important to take the 
limit in which the infrared regulator (the fiducial cell) is 
removed. This step affects the observed difference. In- 
deed, taking states that correspond to the same universe 
but with different fiducial cells, one can see that the ef- 
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feet of the compact remnant gets removed once the cell 
V tends to S t . As a consequence, in the limit when the 
regulator is removed, both the sLQC and the sMMO pre- 
scriptions can be considered to converge in the physical 
sense discussed here (focusing the attention on the kind 
of observables that we have introduced, constructed from 

The existence of nontrivial differences shows that the 
prescriptions are truly physically different, and the dif- 
ference cannot be canceled out by a change of repre- 
sentation. This fact has far going consequences, since, 
contrary to statements commonly found in the literature 
[3~t| , the classical effective description of the system does 
depend on the details of the quantization, and the charac- 
teristic effects of the particular prescription that is used 
have to be taken into account in the process of arriving to 
that description and determining its domain of validity. 
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Appendix A: WDW model 

In this appendix, we describe a geometrodynamical 
analog of the system considered in this paper. This ge- 
ometrodynamical model, built via a WDW quantization, 
has been discussed extensively in the literature (see for 
example Ref. Q). Here, we just summarize the prop- 
erties necessary to define the WDW limit of the LQC 
states. 

The model is constructed following a process similar 



to that of the loop quantization (see Sec. Ill B[ ). The only 
difference is that now the geometry degrees of freedom 
are quantized using a standard Schrbdinger representa- 
tion. The kincmatical Hilbert space is given again by a 
tensor product, Kkin = ^'X'Hkin' where "H? is the space 
defined in Sec. Ill B II and the gravitational Hilbert space 
is now Hkin = L 2 (R,dv). The triad operator still acts 
by multiplication, p\v) = sgn(u)(27T7^p 1 v / A|f |) 2 / 3 |v), but 
the connection is now a well defined derivative operator, 
c = 2i(27T74i) 1/3 A- 1 / 3 |w| 1 / 6 a u |w| 1/6 , contrary to the sit- 
uation found in LQC. Then, the evolution operator ana- 
log to (with a factor ordering compatible with that 
chosen for the latter operator) takes the form: 

<3 = -UnGy^dMdvVWl (Al) 

This operator is essentially selfadjoint in HklT ■ ^ s spec- 
trum is positive, twofold degenerate and absolutely con- 
tinuous. Opposite orientations of the triad (v > and 
v < 0) are disjoint under the action of 0, therefore the 
restriction to the (anti)symmetric sector can be imple- 
mented by considering only the part v > and proceed- 
ing then to the (anti)symmetric completion of that part. 
In the symmetric sector, there exists an orthonormal ba- 
sis of generalized eigenfunctions {e k | of whose elements 
are the rescaled plane waves 

e ±k (v) = (e ±k \v) = 1 e ±jfch - v £ K + . (A2) 

V 27TW 

The corresponding eigenvalues are uj 2 = Yl-nGk 2 . These 
generalized eigenfunctions satisfy the normalization con- 
dition 

(e k \e k ,) = 5(k-k'). (A3) 

The group averaging procedure is straightforward to ap- 
ply in this case, and provides the Hilbert space of physical 
states Kphy = - £ ' 2 ( K > dfc ) 3 where 

*(»,0)= / dk±(k)e k (v)e luj(k)cl> (A4) 



and tu(k) = VUnG\k\. 
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